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ABSTRACT 

A  superconducting  field  magnet  is  designed  and  construc- 
ted for  use  in  a  small  synchronous  generator  of  approximately 
10  horsepower.  The  equations  describing  the  field  flux  dis- 
tribution are  developed  and  compared  with  experimental  results. 
The  magnet  is  tested  to  insure  superconducting  operation  in 
the  generator  system  Is  possible.  The  performance  of  the  magnet 
in  the  production  of  a  superconducting  field  was  found  to 
be  as  predicted.  For  the  construction  of  such  fields  in  the 
future,  an  optimization  scheme  is  presented  for  the  selection 
of  field  geometry,  as  well  as  means  for  predicting  stresses 
in  the  support  structure. 
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LIST  CF  SYMBOLS 


B  Flux  density,  Gauss 

H  Field  Intensity,  ampere- turns /meter 

I  Field  current,  amperes 

J  Armature  current  density,  amperes/cm 
a 

2 

J  Field  current  density,  amperes/cm 

L  Overall  magnet  length,  inches 

Nf  Number  of  turns  in  field  vrinding 

R  Inside  field  radius,  Inches 

R?  Outside  field  radius,  inches 

R.  Inside  armature  radius,  inches 

R  Outside  armature  radius,  inches 
o 


6       Pole  face  angle, 
X      Packing  factor,  conductors/cm 
p  Radial  distance,  inches 

0  Clockwise  angular  measurement  from  magnetic  axis 

of  field  vrinding, 
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INTRODUCTION 

Recent  investigations  into  the  use  of  supercon- 
ductors in  electrical  machines  have  indicated  that 
generators  can  be  constructed  which  will  have  higher 
output  at  greater  efficiency,  lower  weight  per  unit  of 
power  and  lower  synchronous  reactance  than  present  gen- 
erators. Future  requirements  for  high  capacity  turbogen- 
erators will  demand  the  projected  capabilities  of  such 
machines,  and  their  application  in  gas  turbine/electric 
drive  propulsion  plants  for  marine  use  may  successfully 
overcome  the  present  excessive  weight  problem  of  such  plants. 

The  most  promising  use  of  superconductors  in  a  genera- 
tor is  in  the  field  windings,  where  the  current  is  direct. 
Elimination  of  resistive  losses  in  the  field  should  increase 
generator  efficiency  measureably  by  reducing  the  power  for 
excitation.  The  use  of  superconductors  in  the  armature  is 
not  so  promising,  due  to  a  number  of  problems  associated  with 
the  use  of  superconductors  with  alternating  current. 

Since  present  design  philosophy  is  content  to  leave  the 
armature  at  ambient  temperature,  the  remaining  system  design 
consideration  is  whether  the  armature  or  the  field  should 
be  the  rotating  member.  For  large  generators,  the  answer 
lies  in  the  amount  of  power  which  must  be  transferred  to  or 
from  the  rotating  member.  If  the  armature  rotated,  the  amount 
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is  large,  and  the  required  brush  contact  area  is  large. 
However,  if  the  field  rotates,  as  is  the  case  in  present 
large  generators,  the  power  and  brush  area  are  small.  Thus, 
the  most  practical  configuration  is  with  the  field  rotating, 
preferrably  inside  the  armature  to  ease  armature  heating 
problems . 

In  order  to  uncover  impediments  to  the  construction  of 
large  superconducting  machines,  a  small  superconducting 
generator  project  was  undertaken  at  the  Massachusetts  Insti- 
tute of  Technology.  The  project  includes  the  design  and  con- 
struction of  both  the  rotating  dewar  to  house  the  supercon- 
ducting field  at  liquid  helium  temperatures  and  the  field 
magnet,  as  well  as  the  future  design  and  construction  of  a 
shielded  armature.  Figure  1  is  a  schematic  diagram  of  the 
proposed  system.  Presently,  the  dewar  and  field  assembly  are 
intended  for  use  with  the  stator  from  a  two-pole,  three- 
phase  Induction  motor,  rated  at  7  horsepower.  This  thesis 
is  concerned  with  the  design  and  construction  of  the  required 
field  magnet,  subject  to  the  objective  of  producing  the 
maximum  possible  field  intensity* 

Much  of  the  work  done  in  this  thesis  was  based  on  a  pro- 
ject conducted  by  the  Avco-Everett  Research  Laboratory  (1), 
under  which  a  small,  stationary  field  alternator  was  con- 
structed. The  analysis  used  in  that  project  was  extended 
to  the  present  investigation,  which  deals  with  the  inherently 
more  practical  rotating  field  machine. 
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Figure   1.    Proposed  Superconducting  Generator 
System 
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DESIGN  SELECTION 

Since  the  field  iragnet  required,  by  the  design  was 
to  be  incorporated  with  existing  portions  of  the  proposed 
generator  system,  certain  design  dimensions  were  immediately 
fixed.  These  included  the  length  and  outside  diameter  of  the 
magnet.  Additionally,  the  following  demands  were  made  of  the 
designi 

1.  It  was  desired  to  keep  the  superconducting  material 
within  liquid  helium  throughout  the  period  of  operation. 
Since  the  configuration  assumed  by  the  surface  of  the  liquid 
helium  during  rotation  would  be  a  steep  parabola,  all  of  the 
conductor,  including  the  end  turns,  would  have  to  reside  in 
an  annulus  close  to  the  dewar  wall.  y 

2.  The  dimensions  of  the  stator  to  be  used  limited  the 
outside  diameter  and  length  of  the  field.  Because  the  liquid 
helium  would  not  be  replenished  during  the  operation  of  the 
rotor  for  the  Initial  trials,  an  annular  thickness  of  1" 

was  selected  to  take  advantage  of  a  reasonably  long  operating 
time,  and  still  allow  a  signigicant  volume  of  conductor  to 
be  fitted. 

3.  Since  the  predicted  field  intensities  for  supercon- 
ducting magnets  are  so  large,  the  requirement  for  the  use 

of  ferromagnetic  materials  as  flux  paths  has  been  eliminated. 
Therefore,  unlike  conventional  machines,  the  forces  generated 
by  the  field  are  carried  by  the  conductors,  rather  than  by 
any  iron  core,  and  the  structure  had  to  be  configured  to  trans- 
fer the  torque  producing  forces  from  the  conductors  to  the  shaft. 
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kt    In  order  to  keep  the  conductors  at  the  required  temp- 
erature for  superconducting  operation,  adequate  contact  with 
the  liquid  helium  must  be  maintained.  This  implied  that  pass- 
ages would  have  to  be  provided  in  the  structure  to  allow 
liquid  helium  to  reach  the  conductor  from  the  original  intro- 
duction point  on  the  central  axis. 

Consideration  of  several  general  design  alternatives 
led  to  the  configuration  illustrated  in  Figure  2,  The  selec- 
tion of  the  number  of  fins  to  be  used  and  the  dimensions  re- 
quired is  covered  in  Appendix  A.  The  two-pole  field  is 
produced  by  a  superconducting  wire  winding  on  a  cylindrical 
stainless  steel  tube  with  an  outside  diameter  of  3M.  The  wire 
is  wound  into  segments  created  by  fins  extending  radially 
outward  from  the  tube.  Each  segment  contains  approximately 
900  strands  of  .020"  diameter  wire.  The  wire  is  a  .011" 
diameter  core  of  niobium-^0^  titanium  surrounded  by  copper 
and  Formvar  insulated. 

As  shown  in  Figure  2,  the  field  is  wound  in  six  coils, 
each  with  approximately  900  turns.  The  winding  was  accomplished 
in  place,  in  order  to  achieve  the  maximum  space  factor.  Tne 
design  assures  that  all  of  the  conductor,  including  the  end 
turns ,  is  located  in  an  annular  volume  1"  thick.  Circular 
perforations  in  the  central  structural  cylinder  allow  liquid 
helium  to  pass  to  the  conductors  for  cooling.  These  passages 
cover  approximately  20$  of  the  cylinder  surface. 

The  entire  field  is  enclosed  by  a  thin- walled  stainless 
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steel  tube  to  prevent  radial  movement  of  the  conductors  and 
to  assure  that  helium  vapor  produced  by  boil-off  will  be 
confined  to  the  region  inside  the  liquid  annulus . 

The  torque  is  transmitted  by  the  radial  fins  from  the 
conductor  to  the  central  cylinder  and  then  to  the  dewar  and 
shaft  assembly  by  J"  cap  screws.  The  estimated  torques  and 
stresses  experienced  by  the  conductors  and  structure  are 
computed  in  Appendix  C. 

The  anticipated  number  of  ampere- turns  and  the  resulting 
field  intensities  which  could  be  achieved  by  the  design  are 
presented  in  Appendix  B, 

The  surface  of  all  structural  portions  of  the  magnet 
is  covered  with  a  varnish  type  insulation  to  prevent  shorting 
the  winding  bo  the  structure.  Additionally,  the  fins  are 
wrapped  with  a  Teflon  and  fiberglass  tape  to  provide  more 
insulating  protection.  Specifications  for  the  wire  used  in 
the  construction  of  the  magnet  are  given  in  Appendix  D. 

Since  the  final  measure  of  effectiveness  of  the  magnet 
as  a  system  component  is  the  torque  which  it  can -produce,  an 
optimization  scheme  is  presented  in  Appendix  E,  proposing 
that  there  exists  a  particular  geometry  which  will  develop 
the  maximum  torque.  The  results  of  such  a  scheme  will  provide 
a  check  on  the  effectiveness  of  the  thesis  design. 
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Figure  2.  Two-pole,  superconducting  field  magnet. 
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EXPER IMEMTAL  PROCEDURE 

In  order  to  evaluate  the  performance  of  the  magnet,  the 
flux  pattern  was  measured  and  compared  with  that  predicted 
by  the  theoretical  values  developed  in  Appendix  B.  The 
measurement  was  made  at  room  temperature,  because  the  field 
intensity  is  essentially  independent  of  the  environment,  as 
long  as  magnetic  material  is  not  present. 

Since  the  winding  resistance  was  found  to  be  200  ohms, 
the  flux  was  mapped  at  a  low  value  of  field  current  to  pre- 
vent overheating  the  winding.  To  produce  .20  amperes  in  the 
winding,  *J-0  volts  were  required  at  the  terminals.  This  voltage 
was  obtained  from  a  110  volt  d.c.  source,  and  the  excess  power 
was  dissipated  in  wire-wound  and  ceramic  resistors,  as  shown 
in  Figure  3. 

The  pattern  was  mapped  by  placing  the  probe  of  a  rotating- 
coil  flux  meter  in  the  field  and  reading  the  flux  density 
directly.  Readings  were  taken  along  three  loci,  as  in  Figure 
^»  radially  from  a  pole  face,  circumf erentially  at  the  outside 
radius,  and  axially  at  the  pole  face  and  the  outside  radius. 

To  insure  that  the  magnet  would  produce  a  superconducting 
field,  a  static  test  was  made  to  determine  the  resistance  of 
the  winding  in  a  bath  of  liquid  helium.  The  winding  was  mounted 
in  the  dewar  and  energised  with  a  small  field  current  while 
submerged  in  liquid  helium.  The  voltage  and  current  at  the 
field  terminals  were  measured  using  the  setup  illustrated  in 
Figure  5. 
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Figure  3.  Circuit  used  to  map  flux  pattern 


10- 


S>N 


FLUX.  METER 


Figure  k.    Loci  for  mapping  flux  pattern 


-11- 


(i 

D 

^-  ■  ■  I                          

Liquid  < 
Fill 

VACUUH- 

TU8E 

1 

.  .  ,y 

"\  ■ 

VOLTM£T£JZ 

X 

o 
5 

CURRENT 

+ 

Source 

d.o. 

He 

^ 

Level 

5 

S/ 

X/ 

_ 

-|l|i|h 
6v 

+       J 

^ 

* 2 

<^ 

Figure  5.  Apparatus  for  measurement  of  field  resistance 
at  liquid  He  tenperature. 


-12- 


RESULTS 

The  flux  density  map  is  presented  in  Figures  6  through 
8.  The  pattern  along  the  radial  axis  is  shown  in  Figure  6; 
the  measured  values  are  dots  plotted  with  the  predicted 
curve  as  derived  in  Appendix  B.  The  circumferential  pattern 
is  similarly  displayed  in  Figure  7.  The  axial  flux  distribu- 
tion in  Figure  8  is  a  curve  faired  through  the  experimental 
values  of  flux  density. 

The  variation  of  electrical  resistance  with  environ- 
ment temperature  for  the  winding  is  shoi\m  in  Figure  9.  Three 
measurements  were  taken i  one  at  room  temperature,  the  second 
at  liquid  nitrogen  temperature,  80  K,  and  the  last  at  a 
temperature  assumed  to  be  less  than  9-7°K.  The  values  of 
resistance  obtained  are  plotted  against  a  theoretical  resis- 
tance curve  which  is  divided  into  two  parts.  Above  9.7  Kf 
the  theoretical  resistance  curve  assumes  that  copper  is  the 
sole  current  carrier,  and  the  curve  is  based  on  a  resistance 
versus  temperature  curve  for  copper  from  Scurlock  (2).  Below 
9.7°K,  the  transition  temperature  of  Nb-Ti,  the  curve  assumes 
that  the  Nb-Ti  core  is  superconducting  and  carrying  all  the 
current. 

The  results  of  the  optimization  scheme  of  Appendix  E 
indicate  that  for  the  given  number  of  conductors  (5500)  and 
the  given  outside  radius  (2.5" )»  a  pole  face  angle  of  10 
and  an  inside  radius  of  1.875"  would  develop  the  maximum 
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torque  when  used  with  an  armature  whose  current  density 

2      ,   2 

Is  10  amps/cm  . 

The  results  of  the  structural  stress  analysis  In  Appen- 
dix C  show  that  the  scantlings  selected  for  construction  are 
sufficient  to  carry  the  electromagnetic  and  rotational  force 
loads.  The  maximum  longitudinal  bending  stress  anticipated 
in  the  support  tube  is  6200  psi  at  the  mid-length.  The 
maximum  bending  stress  expected  at  the  fin  root  Is  1000  psi. 
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DISCUSSION  OF  RESULTS 

As  the  flux  maps  show,  there  is  essential  agreement 
between  the  theoretical  and  measured  distribution  of  flux 
density  around  the  field  winding.  The  axial  distribution, 
which  was  not  theoretically  derived,  is  of  interest  because 
it  allows  some  insight  into  the  determination  of  the  effec- 
tive length  of  the  winding.  The  shape  of  the  curve  indicates 
that  the  flux  density  is  within  75%   of  the  maximum  value 
over  50%   of  the  length  and  within  90%   of  the  maximum  over 
k0%   of  the  length.  Thus,  the  effective  length  might  well 
be  expressed  as  that  length  over  which  the  flux  density  is 
within  a  given  percentage  of  the  maximum. 

Since  the  resistance  of  the  winding  followed  the 
expected  pattern,  it  gave  encouragement  to  the  ability  to 
achieve  superconducting  performance  with  the  magnet  as 
designed. 

The  fact  that  the  torque  would  be  optimized  for  a 
smaller  value  of  pole  face  angle  and  a  larger  value  of 
Inside  radius  of  the  winding  than  were  possible  under  the 
project  constraints . indicates  that  the  design  of  a  generator 
system  should  Include  the  field  optimization  early  in 
the  design  procedure.  The  existence  of  an  optimum  geometry 
should  be  sufficient  motivation  for  careful  field  design, 
if  the  most  efficient  machine  is  to  be  realized. 
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CONCLUSICMS 

The  field  winding  developed  under  this  project 
performed  satisfactorily  when  examined  for  the  production' 
of  flux  intensity  and  for  superconducting  operation.  The 
geometry  of  the  field  was  not  the  optimum  for  the  number 
of  conductors  used,  since  constraints  of  existing  equipment 
prevented  optimization.  The  design  produced  under  this 
thesis  is  considered  adequate  to  meet  the  requirements 
of  producing  the  maximum  field  intensity  for  the  generator 
system  under  development « 

The  maximum  field  intensity  which  can  be  produced  is 
predicted  to  be  JO   kG,  using  a  field  current  of  80  amperes. 
For  this  field  strength,  the  stresses  experienced  by  the 
structure  will  probably  produce  deformations  well  below  the 
elastic  limit.  The  ease  with  which  superconducting  operation 
was  achieved  is  interpreted  as  a  prediction  that  such 
operation  will  be  possible  in  rotation. 


-20- 


RECOMMEMDATIONS 

Several  experiments  should  follow  naturally  from  the 
work  done  in  this  thesis.  The  first  wou3.d  be  to  increase 
current  to  the  nonrotating  magnet  during  superconducting 
operation  in  order  to  find  the  maximum  field  strength  and 
quenching  current.  Next  the  magnet  should  be  rotated  to 
Insure  that  it  is  properly  balanced,  and  any  balance  errors 
corrected.  Then  a  stator  might  be  fitted  with  a  single  turn 
coil  to  determine  induced  voltage  when  the  field  is  rotated 
in  the  stator.  Finally  the  field  should  be  tested  in  a 
distributed,  three-phase  winding,  in  operational  conditions 
of  loading. 

For  larger  machines  in  the  future,  it  is  felt  that  the 
windings  should  be  made  in  hanks  which  can  be  fitted  later 
into  the  slots,  rather  than  being  wound  in  place.  Although 
there  may  be  a  slight  reduction  in  the  space  factor  which 
can  be  achieved,  the  savings  in  construction  time  and  increased 
simplicity  are  more  than  offsetting. 

When  it  becomes  feasible  to  provide  continuous 
replenishment  of  liquid  helium,  the  inside  radius  of  the 
winding  can  De  reduced,  providing  greater  volume  for  conduc- 
tors and  higher  field  intensity. 


NO  ORDER  SLIPS  INVOLVED 
IN  TECHNICAL  PROCESSING 
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APPENDIX  A 


Selection  of  Design  Dimensions 
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The  dimension  1  f  will  determine  leff ,  the  effective 
length  of  the  field,  so  it  was  felt  that  1  -  should  be 
at  least  as  long  as  the  existing  armature  effective  length, 
or  approximately  3.5" •  Li  the  overall  length,  was  reduced 
from  9"  to  8"  for  the  following  calculations  to  account  for 
the  fact  that  the  wire  could  not  be  wound  to  the  ideal 
shape  by  hand.  The  relationship  between  the  length  required 
for  end  turns  let>  the  overall  length  L  and  the  pole  face 
length  lpf  becomes  let  =  J(L  -  lpf)  =  2.25". 

The  conductor  wires  from  three  slots  must  be  able  to 
pass  through  a  cross  section  1"  x  le+-H  t    so  the  problem  is 
simply  to  select  a  pole  face  angle  6  (Figure  A-2)  which 


satisfies  the  equation 

2 
[2 


7L  (^  -R*)  90   -*  -  let 
IjT   *   J-    90 


Figure  A-2.  Geometry  for  calculation  of  pole  face 
angle. 
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APPENDIX  B 


Calculation  of  Flux  Densities 
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Figure  B-l  illustrates  the  nomenclature  to  be  used 
in  the  development  of  expressions  for  the  flux  density 
distribution.  It  is  first  assumed  that  the  field  current 
distribution  can  be  described  by  current  sheets  of  thickness 
dr  as  shewn  in  Figure  B-l,  These  sheets  have  length  leff % 
the  effective  length,  and  traverse  the  circumference  from 
i    to  (180-i). 

The  current  distribution  can  be  expressed  in  a  Fourier 
series  determined  as  follows » 

K(r)  =  21  kjj  sin  mx&  L  a  7\ 

n  odd       L 
,x    . 

k  =2   K(r)  sin  n^zfd^ 

11    7T  Jo 


•  18C-i 

^n  =  —  \  Jf  d"r  sin  n^^ 


K(r)  =  XI      2_Jf  dr  [-cos  n(180-S)  +  cos  nSJ  sin  n0 
n  odd  n  t 

Once  the  current  distribution  is  known,  the  field  in- 
tensity may  be  found  by  solving  Laplace *s  Equation  in 
spherical  coordinates, 


2 
?   dp    v  dp'  r     (30 

where  the  field  intensity  H  =  VU,  such  that 


1   A-.  ((sM)      +     lz     *    U2      =   ° 


H^  =  J_U  and  H0  =   1    3U 

The  solution   is   subject  to   the  following  boundary   conditions i 

1.  H  is   finite  at      ^  =   0 

2.  H  ->    0  as    p  ->  &o 
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Figure  B-l.  Peripheral  distribution  of  surface  current 
density  in  sheet  of  thickness  dr  for  two- 
pole  field  winding. 
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3.  At  the  current  sheet,  H*  is  continuous,  i.e., 

H^  (^>r)  -  H^  (p<r)  -  0 
k.   At  the  current  sheet 

H^  (/o>r)  -  Bp   (^<r)  =  K(r) 
A  product  solution  is  assumed  with  the  formi 

/       n  ~nx / 

U  =  (A  jo  +  Bo     )(C  sin  np   +  D  cos  n0  ) 

Applying  the  boundar;/  conditions,  the  solution  becomes  i 

f  or  fo  <  r ,   U  =  (  r~*    x>  )  ( -  21  2  Jf  dr  cos  6   )  cos  n  & 

n  n   odd   n  it 

for  x3>r,      U  =    (   r         p     )  (-  21     2  Jf  dr   cos  %  )cos   n0 

n         '  n  odd  ni( 

Since  H  =    VU,    the   increment   of   intensity  at  P  is 

for^o>r,    dH     =    21   _2Jf  dr  r  *    '    P         cos  6  cos   n0 

n  odd"    nij 

for^>>rs    dHw  =    21   2Jf   dr  r]1+1   o"n"1cos  ^  sin  n0 

n   odd-n  ^ 

t-    r,  -n+1     n-1         i 

for  /o<r,   dH*  =    2-    2Jf  dr  r  o        cos  d  cos   np 

n  odd-  n  ^r — 
for/o<r,   dHv  =    21  2  Jf  dr  r"n+1  />n~1cos  6  sin  njzJ 

n  odd-^-         r 

The  field  intensities  are  then  found  by  integration  over 
the  current  carrying  annulus :  ^ 

for  o>R0,  H/>  =  21  2J  /o"n"   cos$  cos  n^  \  rn+  dr 

1         *        r     n   odd-rrV-  '*. 

H0  =     21  2J    /o"r        cos<$  sin  n£    \    rn""  dr 

n  odd-H^r  ^ 

for  Ri</o<R2,  •  *t 

H     .      21  2Jf      cos  t  cos   nJp31"1  )    rn+1dr  +^"1  j    r"n+1di 
n  odd-^Y  r         ^  \ 

_  0  k  f  -n-l  f  *n+l  n-1   f     -n+1 

H     -       ^-     Jf_  cos  o  sin  n^    p         \    r       dr   +p  r  di 

n  edd  r  7^  L  '  ^  /  /^ 
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for  /o<R1,    H/0  =    £  g-Jf/o""1   cosi   cos   n0    (  r"R+1dr 

n   odd  n  7r  y^ 

njzs   -    2w   ZJ.r/Q         cos  6  sin  nrf   \  r  dr 

n   odd   nir  ) 

After   integration,    uhe   equations   become j 

/o<R1»    R.   =]T     2 Jf  ^n"      (R~n+  -   R  )    cos  j  cos   n? 

n   odd   nV(2-n)      *  1 

H0   =  23     ZJ^./o""1    (R:n*2-  R7n+2)    cosi  sin  n^ 
n   odd  nxTf  (2-n)      *  1 

R1<o<R2? 

v-n+2  .     vn+2 


H,o  =    H     2Jf/°         2    f-2n   +    (2+n)/Wj  -    (2-n)/R^       JcosScos 
n  odd  nn  (4-n    )L                         V  /  If1/       J 

Hp  -   72    2J./Q        ,2    [4  +   (2+n)(RVr  -    (2-n)/RAn+2lcos^sin  ] 
n   odd  n^  (4-ri    )                          l^2/  (p1)       J 


v— *  n    1        n+2       n+2 

R2</°  '      H,o  =    2^   2J*/>"  <R?      -  Hi       )    cos  5  cos   np* 

'  r     n   odd  nV(n+2)  * 

H0  =    2_J    2Jfv^  (Ro     -  Hi      )    cos&   sin  n0 

n  odd  nVuT+2)         X 

To  evaluate  these  equations  for  the  design  selected,  Jf 

must  be  computed.  Since  the  total  ampere-turns  are  equal  to 

the  current  density  integrated  over  the  sectional  area  of  the 

conductor  volume,  the  ampere-turns  are 

2    2 

Nf  if  =  teJL&zjzJLil  12°^<&.     ' 

1    1  Z  90° 

For  the  design  as  constructed,  N~  =  5500 ,  R^  =  3.82  cm, 

t  ° 

R«  =  6.35  cm,  b  -   25.5  t  These  values  may  be  used  to  solve 

2  2 

for  J~ ,  which  becomes  J-  =  1.98  10  I  amp/ cm  .  It  is  most 

informative  to  obtain  the  ratio  of  flux  density  to  field  current 
for  the  field,  so  the  relationship  B  =/^0H  is  employed.  With 
this  relationship,  the  value  of  Jf   computed  above,  and  the 


-30- 

field  dimensions,  the  field  int?nsity  equations  may  be  trans- 
formed to  a  description  of  the  flux  density  distribution  for 
the  field  magnet.  Since  the  measured  flux  distribution  is  to 
be  compared  with  these  distributions,  it  is  sufficient  to  con- 
sider only  the  fundamental  component  of  the  flux  (i.e.  n  =  1). 

Thus,  for  /o<R1»  B^  =  2  Jf //„  (R2  -  R1)  cos  S  cos0 

fr 

B0   =  2J^>  £<«,  (R2  -  R^)    cos  b  sin^ 

7T 

3 

for  R1<)o<R2,    B^   =   2Jr//.o[:2    +  ^2   "(^l)  J  cos  ^  cos^ 

B0  =  2 Jf/y ft  -  2R2  -(Ri)  ]  cos&sin^ 

3         3  £ 

for  R2</0    i      Bf  =  2Jfi^2    (B2   ~  Rl  ^    C  C0S^ 

3f/> 

B<*  =  2Jf£0  (R9  "  R^  cos^  sin^ 

If  these  equations  are  evaluated  at  the  pole  face  (0  -  0), 
an  expression  for  the  radial  flux  distribution  is  obtained: 
f or  p<Rj   B^  =  3.60  102  If  Gauss 

for  R1<p<R2  Bp  =  ^8.2  If  [-2p+  19.1  -  ii*I2]  Gauss 
for  R2<p   B^  9.60  |03  If   Gauss 

where  p   is  in  centimeters.  The  circumferential  distribution 
at  any  radius  can  be  obtained  by  using  the  above  values  as 
maxima  and  applying  the  relationship  BA$)   =  B/>max  cos^. 
The  values  so  obtained  will  be  a  reasonable  approximation 
over  the  effective  length  in  the  axial  direction. 
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APPENDIX  C 


Calculation  of  Principal  Stresses 


in  the  Field  Assembly 
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The  stresses  In  the  field  assembly  can  be  convieniently 
separated  into  two  catagories,  those  arising  from  the 
electromagnetic  forces  and  those  due  to  the  rotation,  and 
then  superimposed  to  yield  the  total  stress.  Of  greatest 
interest  are  the  stresses  in  the  central  support  tube  and 
in  the  cantilevered  fins.  The  first  examination,  then,  is 
the  determination  of  the  stresses  in  the  central  support 
tube  due  to  the  electromagnetic  forces.  These  forces 
originate  in  the  conductors  and  are  transmitted  to  the  struc- 
ture, much  like  a  pressure  distribution  over  the  surface. 

The  force  on  a  differential  volume  in  the  radial 

direction  is 

dFJo  =  J  x  B,dV 

With  the  expression  for  the  fundamental  component  of 

tangential  flux  density  derived  from  Appendix  B  and  (1), 

the  differential  force  becomes 

2  3  3     3  1 

a^gfggJf  [-2+2R2-/R^]|0sin^cos£   -J^Jf  Ja(R0-Rj.)  (R2-Ri)cos*jdV 

The  force  distribution  per  unit  length  as   a  function   of    0 
may  be  found  by   integrating  the  above  differential  over   the 
thickness    of    the   conductor  volume.    The  result  of   this    inte- 
gration is 

Ff  (0)    =  lMJt  [-2(R2-R1)    +  2WV   "  RjlnfRg-R^JcosSsIrij* 
-  gk  Va    (VRi)(R2"Rl}    C0S^   ' 
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2    .  2 

If  a  value  of  10  amp/cm  is  assumed  for  armature  current 

2  2 

density  and  1,98  10  I  amp/cm  for  field  current  density 

as  found  in  Appendix  B,  then  for  the  geometry  of  the  gen- 
erator system  anticipated  by  the  MIT  project,  where 
R  -  Ri  =  2,5^  cm,  the  force  distribution  on  the  outer 

surface  of  the  central  support  tube  is  given  by 

2 

F^  =  .0617  If  sin£  -  ,0396  If  cos0 

From  Appendix  B,  the  maximum  field  density  is  given  by 

2 

B/>  =  3r60  10  I  .  If  this  relationship  is  plotted  on  the 

short  sample  performance  characteristic  of  Figure  D-l,  the 

intersection  gives  the  maximum  allowable  field  current, 

approximately  80  amperes.  With  this  value  of  Ifl  the 

force  can  be  closely  approximated  by  the  first  term  only, 

Fx,(0)  =  39^  sin0    for  O<0<TT 

-39*4-  sin?        for  if<  0<  2  tt  ,  in  newt  ons /radian. 

From  Wang  (3)f  the  equilibrium  equations  for  circular 

cylindrical  shells  are 

2        2  2  2 

d  u  +  1-g  3 u?  +  l+v  _<3_v  -  v_  3  w  +  p.  (1-V   )   =  0 
a?   2a2"  a 0"  2a  d*30  a8>"       Eh 

2         2        2  2    3     3  N 

1+z;  a  u  +  1^  d  v  +  l^itx.  -  I?  5ii  +  h  ?(  d„w  +  ^  ttq ) 

2a  dzty        2   £>x2   a^302   o.     d0      12a  Vdx^  a^07 

+  Jl.JL±  Ax~  +  i2  Azo")  +  Poii^l    =  ° 

2&*\   2      dxz       a       djzJV 


12a^\    2       ar  a       djzJV  c      Eh 

2        3  £                    ^  *     \ 

y_^u  +    3v     -  w  -  h  r  fa-^  c?  W),  +  2a  c)  w  +      c)  w,,  ) 

Ay        p^Gf          a      12a^\       r)  x             d~P-bd>  a  d<25  v 


2        3  ^  zj.  4 

1  r  f  a^  c?  w,    +  2a  3  w       +      3 
la^V     "Jx  cbc2^       ac 

3  3 

h_  /  1   c)  v     +  U    S_v«0\      +  aP<JL 
12\aax2d0       aJ    a^-^y  J     Eh 


-3-  -  'i-vh     =  0 
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where  u,  v,  and  w  are  deflections  in  the  x  (axial),  0 
(tangential),  and  p   (radial)  directions,  respectively. 
If  the  ends  of  the  cylinder  are  assumed  to  be  simply 
supported,  then  the  boundary  conditions  at  the  ends  may 
be  expressed  by  v  =  0,  w  =  0,  and  the  shear  force  and 
bending  moments  are  zero  at  x  =  0  and  x  =  L.  Thus  a 
solution  of  the  form  of  three  Fourier  series  may  be  used 
to  describe  the  deflections i 


u  as  /  .  /      A   cos  n0  cos  mjrx 

m=l  n=0     mn  L 


v  =    52  Yli    B       sin  rift  sin  mttx 

m=l  n=0     mn  L 

w  =   12  22     cmn   cosn^  sin  nmy 

m=l  n=0  L 


The   loading  may  also  be  described  by  a  Fourier  series    obtained 
in  the  usual  wayi 

Pt    =   0,    p~  =  0,    p.   =   394  sinjzf  for   0<#<tl 

1                c  J     -394  sin$rf  for  ir<<2<2% 


oo     «** 


Po  =  S  2C  Dmn  cos  ¥sin  HIS 
J    m=l  n=0  L 

with    D^  =  8    \  394  sin/^cos  nf  dp 

Evaluation  of  the  integral  leads  to  the  series 

d  =  X^  T*  3152   cos  nd  +  1    cos  ncf  sin  m£x 

■*    m  n  mi^-   (l-n)(l+n)  l* 

Substitution  of  the  series  expressions  for  u,  v,  w,  and  p^ 
into  the  equilibrium  equations  gives  three  homogeneous 
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equatlons  In  Amn,  B   ,  and  Cmn.  Values  for  u,  v,  and  w  may- 
be obtained  with  sufficient  accuracy  by  evaluation  of  the 
coefficients  at  m  =  1,2  and  n  =  0,1,2;  since  the  series 
converge  rapidly.  The  primary  stress  in  the  central  structural 
tube  will  be  due  to  the  longitudinal  bending  moment,  which 
is  given  by  Wang  (  3)  as 

2  2 

(MY)    =  Eh3    fa  w  +  _v  /ay_  +  3_w_Y| 
max  12Tl-^2  )  L  a  x2    a*2  V  B^    a^J 

This  expression  may  be  evaluated  for  the  particular  appli- 
cation of  this  design,  using  the  deflection  series  values 
and  other  parameters  as  follows: 

E  (modulus  of  elasticity)  =  30  x  106  psi 
h  (tube  wall  thickness)   =  .120  inches 
V    (Poisson's  ratio)       =  .25 
a  (tube  radius)  =1.5  inches 

L  (tube  length)  =  9  inches  , 

resulting  in  a  bending  moment  of  15  inch-pounds.  Wang  shows 
that  the  stress  associated  with  the  maximum  longitudinal 
bending  moment  is 

<r-  =  6  M       psi 

x   h2  -^max 


Thus,  <r    =     6200  psi  at  mid  length, 
x 
The  longitudinal  stress  due  to  the  rotation  of  a  tubular 

shaft  may  now  be  added  to  that  due  to  the  electromagnetic 

stresses.  The  stress  may  be  calculated  using  the  formula 

<T      =  (3-2fc0  Uf>::2      (b2  +  a   -  2r  )  , 

x       W^v)  3-2i/ 
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whlch  assumes  that  the  ends  of  the  tube  are  constrained 
longitudinally.  For  the  present  design  p=    .278  lb/in^, 
6J  =  3600  RPM,  b  (the  inside  tube  radius)  =  I.38  in,  and  the 
stress  is  determined  for  r  =  b.  This  stress  is  found  to  be 
small,  less  than  100  psi,  in  comparison  to  almost  10,000  psi 
due  to  the  electromagnetic  forces.  The  rotational  stress  is 
therefore  neglected,  and  the  -working  stress  of  the  cylinder 
is  estimated  to  be  an  order  of  magnitude  lower  than  the 
yield  stress  of  the  stainless  steel  used  in  the  support 
construction. 

To  estimate  the  bending  stress  in  the  radial  fins,  there 
are  two  sources  of  loading  which  must  be  considered.  The 
first  results  from  the  forces  which  produce  the  machine  torque 
and  the  second  are  the  field  self-forces.  In  analyzing  the 
first,  it  is  assumed  that  the  machine  torque  is  produced  by 
a  uniformly  distributed  load  on  the  fins  (or  an  equivalent 
point  load  at  the  mid-span  of  each  fin).  From  Appendix  E, 

the  torque  is  found  to  be 

3    3 
T  =  2  m  Jf JQ  (E0  -  Ri)(R2  "  Ri)cos^  1eff 

For  the  following  values  j 

2  7 

Jf  =  1.98  10  If  amps/cur 

J  =  10"  amps /cm 
a 

leff  (effective  length)  =  4  in  =  .1  meter 
If  =  80  amps , 
the  torque  becomes  T  =  ^0  ft- lb. 

Since  this  torque  will  be  carried  by  12  of  the  1^  fins, 
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3.3  ft-lb  is  the  torque  on  each  fin.  The  radial  center  of 
each  fin  is  at  a  2  inch  radius,  so  that  the  equivalent  load 
is  a  concentrated  force  of  20  pounds  at  the  inid-radlus  of 
the  conductor  annulus .  The  bending  moment  resulting  from 
this  force  is  .84  ft-lb,  The  worst  stress  condition  will 
occur  at  the  intersection  of  the  shortest  fin  and  the  support 
tube,  since  the  shortest  fin  will  have  the  smallest  moment 

of  inertia.  Figure  C-l  shows  the  geometry  of  the  point 

3 
being  examined.  With  e-=  Mjr  and  I  =  bh  ,  the  stress  is 

I         12 

found  to  be  100  psi,  well  within  the  tolerances  expected 

of  the  weld  and  fin  material. 
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Figure  C-l.  Geometry  at  Fin  Root 
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In  examining  the  fin  loading  due  to  the  field  self- 
forces,  the  differential  of  force  on  the  conductors  can  be 
written 

&Fp    =   J  x  B  dV 
The  force  on  the  conductors  in  the  segment  from  0-  to  0 
can  be  found  by  evaluation  of  the  integral 

F0  =  j  j   \  Jf  *  Bp(odpd*dl  , 
which  becomes      '  ^' 

Fp=\\   \    2^°Jf  -2  +  3£2  "(-l)   (°  cos  ^  cos?    dfd^^l 
o  ^  ni, 
when  ths  expression  for  B  is  inserted.  With  the  proper 

values  of  field  variables  applied,  the  expression  becomes 

2     1^2 
F  =  .951  I  sin^j     newtons . 

This  expression  is  maximized  for  the  fins  closest  to  the  pole 

face,  where  0^  =  25.5  and02  =  ^7  •  The  maximum  value  is 

F  =625  newtons  or  275  pounds.  If  this  force  is  assumed  to 

be  uniformly  distributed  over  the  fin,  the  resulting  bending 

stress  at  the  fin  root  is  1000  psi,  an  order  of  magnitude 

larger  than  the  stress  resulting  from  the  torque  load,  but 

still  well  within  the  expected  limits  of  the  structural 

material. 
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APPENDIX  D 


Superconductor  Specifications 


-*J0- 


The  following  data  are  taken  from  an  advance  tech- 
nical bulletin  describing  Nb-Ti  Supergenic  Wire,  Type 
SG-100D,  produced  by  Avco-Everett  Research  Laboratory. 

a.  Material  -  Supergenic  Wire  is  a  composite  super- 
conducting wire  composed  of  copper  clad  Nb-^0$  Ti . 
This  material  permits  the  economical  manufacture  of 
superconducting  coils  up  to  100  kG. 

b.  Dimensions  and  Weight  -  Type  SG-100D  has  a  0.020" 
outside  diameter  and  contains  a  0.011"  core  of  Nb- 
k0%   Ti,  The  insulation  Is  a  one  mil  coating  of 
Formvar . 

c.  Design  Stress  -  A  design  stress  between  10,000  psi 
and  15 f 000  psi  is  currently  recommended  on  the  total 
cross  section. 

d.  Stability  -  Type  SG-100D  is  fully  stable  in  short- 
sample  tests  and  can  be  incorporated  into  fully  stable 
coils  when  adequate  cooling  is  provided. 

e.  Handling  -  Supergenic  Wire  is  a  rugged,  ductile  material. 
Both  copper  and  niobium  titanium  are  ductile  materials 
and  may  be  handled  with  conventional  winding  machines. 
Lapped  soldered  joints  are  permitted  and  may  be  wound 
into  a  properly  designed  coil. 

f .  Performance  -  Developmental  type  SG-100D  short-sample 
performance  is  shown  in  Figure  D-l.  A  manufacturing 
tolerance  of  -10$  currently  is  typical  of  developmental 
quantities . 
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Figure  D-l 


Supergenic  Wire  Type  SG-100D  short  sample 
performance  characteristic  with  magnet 
operating  characteristic  superimposed. 
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APPENDIX  E 


Optimization  of  Field  Geometry 
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In  approaching  the  design  of  a  superconducting  field 
in  which  none  of  the  variables  of  size  are  initially  fixed, 
it  is  evident  that  higher  field  intensities  will  be  obtained 
as  the  machine  size  is  increased,  since  there  is  no  longer 
any  limitation  imposed  by  saturation  of  ferromagnetic 
material.  The  quenching  field  strength  of  the  conductor  can 
also  be  controlled  by  varying  the  geometry  of  the  conductor. 
Therefore,  it  is  assumed  that  in  practical  machine  design 
(4  ),  the  objective  which  must  be  met  by  the  rotor  is  that 
it  produce  the  maximum  field,  subject  to  a  given  outside 
diameter  and  fixed  volume  of  conductor.  In  the  generator 
system,  this  objective  translates  to  the  requirement  that 
the  torque  produced  be  a  maximum.  Figure  E-l  illustrates  the 
field  geometry  to  be  optimized,  namely  the  values  of  &  and 
R^  are  sought. 

The  torque  can  be  expressed  by  the  integral 
T  =  J  ( J  x  B  x  p)f  df dldl  , 
where  E  is  the  field  produced  in  the  armature  and  J  is  the 
current  density  of  the  field.  From  (1)  the  fundamental 
component  of  the  radial  flux  density  inside  the  armature, 
per  phase,  is  given  by 

B/°  =  3M*   (Ro  "V  cos  t^-^). 
2  7f 

Then     J  x  B  =  3/t^JfJa    (R0  "   Ri^    cos    ^  ~^  ^    in   the  axial 

— -^re- 
direction. 


a 


a 


Figure  E-l.  Nomenclature  for  Field  Geometry  Optimization. 
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Since   Jf   =  +  Jf        for  o   <  0  <  180   -  h 

-   J         for   180   +  &<#<360   -6 

0  elsewhere, 

the  Integral  becomes 

*z  180-8  A^oA 


\  .180-4  f/ 


^1   ©  ^1  I80+-S 


assuming  that  the  maximum  torque  is  delivered  for  f -   90  . 

Evaluation  of  the  integral  yields 

3    3 
T  =  2jU  Jf  Ja  (R0  -  R^Ug  -  Ra  )  cos* 

This  expression  for  the  torque  is  to  be  maximized,  subject 
to  the  following  constraints i 

1.  a  fixed  outside  rotor  radius,  R  =  r 

c  o 

2.  a  fixed  volume  of  conductor,  Nf ,  where 

2  2 

Nf   =    T  A     /90   -  6    URp   -  R-)    ,    with  a  packing  factor  A 

1     2acl   90   / 

and  conductor  cross-sectional  area  a  . 

c 

The  second  constraint  can  be  rewritten  in  the  form 
£  =  7r   -   N^ac  .    2 L2_ 

T"     XA    TR2  -  RJ 

If  kj  =  If/Jf  and  k2  =  l^Ig^  ,  then 

T  (R1,R2,M  =  k2  If  (r|  -  R^)  cos  h   . 
The  constraints  may  be  written  in  functional  notation  as 

where  k_  =  Kfa_  ,  and 
^     A° 

2.  G  (Rlf  R2,  6)  =  R2  -  rQ  =  0 
At  an  extreme  point  of  the  function  T, 
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Tn  dR1  +  T_  dR0  +  Tjrdi   =  0  , 

FR  dR1  +  FR  dR2  +  Fsdi   =  0  ,  and 
GR  dRl  *  GR  dR2  +  G<d^   =  0  . 

The  subscripts  denote  partial  derivatives  with  respect  to 
the  variable  indicated.  If  the  required  operations  are  per- 
formed on  the  functions  T,  P,  and  G,  and  the  resulting 
values  substituted  into  the  partial  derivative  equations ,'a 
set  of  homogeneous  equations  is  obtained  which  may  be  solved 
simultaneously.  The  resulting  equation  expresses  the 

requirement  for  maximum  torque » 

3    3  2    2 

-2k  (R2  -  R^)sin8   +  3%  U2  -  R1)cos6  =  0. 

When  this  equation  is  combined  with  the  constraints  Rp  =  rQ 

and  8  =  7T   -  k~  _o_i_2   »  tne  requirement  is  expressed 
2     J   "(R2  -  RJT 

in  terms  of  the  single  variable  R.j 
-2k3(r^  -  RjScosf  k3  \    +  3Ri(R2  -  R^sln/  ^3  J  =  0  . 

Evaluation  of  this  expression  for  the  outside  rotor  radius 
value  of  2,5"  and  N  =  5500  turns  leads  to  an  optimum  torque 


when  S   =  10°  and  R±   >  1*875  ■ 
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